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ABSTRACT r 


Algebraic method techniques are applied to design linear control 
systems. Computer programs to accomplish the design procedure are 
presented. Several combinations of feedback and cascade compensated 
systems using different specifications are studied. Specifications | 
considered include the system damping ratio, undamped natural 
frequency, bandwidth, steady state error coefficients and root 
sensitivity. 
The method is shown to be applicable to any transfer function and 


capable of working with any number of parameters. 
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I. INTRODUCTION 


The design problem can be stated as "the problem of determining 
a practical set of physical parameters for a system such that a desired 
performance is obtained." [1] 

Many techniques have been devised to accomplish the design of 
control systems: the Root Locus and the conventional frequency response 
techniques being the most popular. All of these techniques present the 
advantage that the designer has a graphical representation of the 
problem and therefore can visualize the effects that a change in one 
parameter produces in the others. However, they are time consuming 
and have the disadvantage of being unable to handle more than one 
parameter. One way of overcoming this difficulty is by using another 
graphical technique capable of handling more than one parameter. For 
this reason Siljak [2] has developed the "parameter plane" method, in 
which two parameters may be varied at the same time and the solution 
may be obtained by using the transformation proposed earlier by 
Mitrovic [3]. 

The theory and methodology developed for the parameter plane have 
provided the feasibility of using a pure analytic multiparameter pro- 
cedure called "algebraic methods," by means of which the system 
differential equations and the specifications expressed in analytical 
form can be reduced to a set of algebraic equations. The resultant 


algebraic equations, however, are nonlinear and in all but the most 


simple cases solution by hand calculations may be a more difficult 


problem than the actual problem of design. This difficulty has been 
greatly overcome in the last decade by the rapid development of digital 
computers, which have made possible the application of numerical 
techniques for the solution of nonlinear equations. 

In this paper two methods are used to solve the design problem. 
The first method uses a linearization procedure to solve the resultant 
nonlinear equations, while with the second method the equations are 


solved by an iterative method. 


II. ALGEBRAIC METHODS 


The basis for Mitrovic's approach to the root finding problem [1] 
was the Cauchy principle of argument which states that "the conformal 
mapping of a closed contour on the S plane through a mapping function 
produces a contour on the polar plane, and the encirclements of the 
origin of the polar plane by the mapped curve can be interpreted in terms 
of the singularities of the mapping function that are enclosed by the S 
plane contour." Mitrovic chose, for the values of S in the second 
quadrant, the transformation of variable: 

S= W exp li( TT /2)+8) = -W, Sin @ + jW_ Cos 8 (2-1a) 
for a mapping function given by the polynomial on the S plane 


n 


ms ame 
F(s) = >. A S* = 0 Co) 
k=o 


If in equation 2-la Wn is the radial distance from the origin to a 
particular root (later on this will be called the dominant root) of F(s), 
and 8 the angle formed between the second quadrant radial line and the 
negative real axis, and also defining cos 8 as the damping ratio | ; 
then equation 2-la can be written as 


_ 2.1/2 
Sw ot iW AE = e) (2-1) 


Furthermore if 2-1 is substituted in 2-2 it can be shown [2] that the 


powers of S formed obey the following relations 


sk =wk a @)) + sta - §7)1/7u, © #9} 


where 


a k 
iG 4) = (-1) TI @) 


UL $) = (-1)*u_. (¢) 


k 


and U, and T are Chebyshev functions of the first and second kind 


k 


respectively and are given by the following recurring relations [2] 


fot (rSae sd ara TF 
hey W326 yt Y) 


with 


T ( $)=U,( ¥)=1 
fe) 1 


U,( ¥) = 0 and T,( $) = ¥ 


fo) 
Using these transformations and equating to zero independently real and 
imaginary parts, equation 2-2 can be written 


n 


S 7 -nawhuy CF) =0 (2-3) 
k=o 

= kK. ..k —" 

i, (-1)"A,W> U,( $) = 0 (2-4) 
k=o 


Equations 2-3 and 2-4 are the basic equations developed by Siljak [2] 

to be used in the "parameter plane" method of design. They are also the 
basis, as will be apparent in the following sections, by which a pure 
analytical design is made. That is, design using pure analytical 
procedures consists basically of solving equations 2-3 and 2-4 for the 
unknown parameters of the characteristic equation formed by the system 


differential equations, reduced to a transfer function in terms of the 


10 


Laplace transform (S domain), and the adequate cascade and/or feedback 
compensators required to satisfy the prescribing specifications. Itis 
obvious that the coefficients AL. in equations 2~3 and 2-4 will bea 
function of the unknown parameters of the system and those introduced to 


it by the compensator devices. If these parameters are called r, 


(i =1, 2, .........-..-Ff) then the coefficients can be written as 


k ki i=1,....,P (2-5) 
Requirement for solution is that the number of equations must equal the 
number of unknowns. Therefore a set of P equations are required to 
attempt a solution. The trivial way of finding these equations is by 
assuming P roots of the characteristic equation of the compensated 
system and replace each pair of complex ones in 2-3 and 2-4, and each 
real one in the equation 


n 


k 
ss AL Oo -=0 (2~6) 
k=o 


Another philosophy consists of substituting in the basic equations 
only the roots whose real parts have the least negative value and 
dominate the time response of the system (by definition this is the 
dominant root) and derive additional algebraic relationships from the 
system performance specifications. In previous work algebraic relation- 
ships for specifications such as bandwidth, sensitivity and error co- 
efficients have been found. However, they have the disadvantage that 
depending on the complexity of the system they lead to a set of algebraic 
equations whose nonlinearity is such that the solution is extremely 


difficult and assistance of a digital computer is mandatory. 


ll 


Ill. THREE SPECIFICATIONS DESIGN 


A. ROOTS AND ERROR COEFFICIENT SPECIFICATIONS 

Three common performance specifications are the dominant roots 
(actually two specifications) and the error coefficient. The cominant 
roots are a measure of the time response of the system while the error 
coefficient can be interpreted in terms of the steady state errors 
(classical interpretation) or as a measure of the moments of the impulse 
response of the closed loop system function [4]. From this last inter- 
pretation, time delay and rise time can also be derived. 

The design problem in terms of the three specifications, or in 
general in terms of any number of specifications, can be broken into 
three parts: 

a. The specifications are interpreted in terms of constraints 
on the mathematical functions describing system characteristics. 

bis A number of systems are determined, each of which satisfies 
the prescribing specifications; and 

oP The "best" system is selected. 

If conventional frequency response and root locus methods are used the 
procedure that the designer follows is to modify and build up the open 
loop transfer function until the closed loop meets the prescribed speci- 
fications (basically trial and error); but if pure analytical procedures 


are used the problem can be solved in a rather fast and accurate way. 
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1. General Statement of the Problem 


Design a linear control system whose uncompensated open 
loop transfer function, in terms of Laplace transform, is given by a ratio 
of two polynomials and whose closed loop must be such that a given set 
of specifications for root locations and error constant coefficient are 
satisfied. The prescribed specifications imply that basically three 
equations with three unknowns must be found in order to attempt a 
solution; that is, the root specifications lead to a set of two equations 
and the other equation is derived from the error coefficient specification. 
The designer therefore must find appropriate compensator devices which, 
when combined with the system equations in a closed form, meet the 
specifications. Several combinations of feedback and cascade compen- 
sators may be used. However, inthis section, a combination of single 
filter cascade and velocity feedback (tachometer feedback) is used to 
develop a computer program capable of handling a rather general problem 
in which the transfer function is assumed formed by the ratio of two 
polynomials expressed either in its factored form or in powers of the 
complex variable s. Since in most devices used for control systems, 
the order of the numerator polynomial is much less than the order of the 
denominator, it is further assumed that the uncompensated denominator 
transfer function is at least two orders greater than the uncompensated 
numerator transfer function. The system block diagram to be analyzed 


can be arranged as shown in Figure 3-l. 


L3 


Figure 3-1 


Cascade and Feedback Compensation for Roots 
and Error Coefficient Specifications Design 


Considering G, (s) = 1, the system transfer function can be written as: 


2 


m m-l 
G(s) =K N(s) _ x A,S +A,8 Fees tA (3-1) 
f D(s) i i 4 


n n= 
Q,5 +Q,5 Fane AQT 


Likewise the cascade compensator (single filter) can be written as: 


+ Z 


G,,'s) . S+P (3-2) 
and the feedback compensator as: 
H(s) = Kt.sS (3-3) 


Equations 3-1, 3-2 and 3-3 contain four unknowns, namely Kf, Z, P 
and Kt, which when combined to form the characteristic equation of 
the system will be a function of the coefficients, i.e., 
A, = G, (KE, Z, P, Kt) (3 -4) : 
In fact if the characteristic equation is defined as F(s), then 
F(s) =1+Go (8) =0 


F(s) =1+ Ge(S) .K, (3-5) 
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G_.(s)G(s) 
cl 


Q 
_— 
iJ) 
— 
I 


: 1+1G .(s)G(s))K S$ 
cl t 


@ 
—_ 
13>] 
— 
i 


K G (s) 
ee 
Defining the closed loop transfer function as C(s\)/R(s), this can be 


written as: 


C(s) _ G (s) _ K G (s) 
= (eo) = e e 
R(s) 126 @) b4 eG (6) 
(eo) ee 


KK N (s) (S+Z) 


Cis) _ is 

R(s) (S+P)D(s) + K kK SN(s) (S+2) + KK N(s) (S+Z) 
C(s) _ Pls) : 
R(s) ~ F(s) —— 


Expansion of equations 3-5 and 3-6 was needed in order to write the 
computer program in Appendix D. However, the author considers that 
no insight, beyond the fact that the coefficients of P(s) and F(s) are 
functions of the unknown parameters, is gained by writing again such 


expansion. 


If n is the order of the uncompensated denominator polynomial 
and m the order of the corresponding numerator and if n= m+2 then the 


order of the characteristic equation will be n+1; therefore, it can be 


written as: 


F@) = G Kk K ZFS" = 0 (3-7) 
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From the basic eauations 2-3 and 2-4 and the given root specifications 


expressed in polar form (equation 2-1) the following equations can be 


written: 
nt+1 k . 
(-1) G, (KK, .2,P)W Uy ¥$ ) = 0 (3~8) 
k=o 
n+l k k 
yy (-1)"G, K,,K,.Z,P)W "UL ( $)= 0 (3-9) 
k=o 


The error coefficient specification can be expressed in algebraic form 
in the following way: 
let Goc(s) be the open loop transfer function of the compensated system, 


then » 


if K-57 + (a ZA Mas 5 one AMZ) 


G (s)= oi a RS Bt ee se ee 
oc 


nt+1 n m+2 e 
(Q,s (Q,P+Q,)S +. tQ) PKK (A,S Hine ame 


From this expression the three more common types of systems can be 
defined in terms of the coefficients in the following way: 
Type zero system: Qt is not equal to zero 


Type one system: Qn is not zero and Qe is equal to zero 
Type two system: or and oF el are both equal to zero 


The constant error coefficient for a type zero system can then be written as 


K ZK ; 
K =limG _(s) = “putts (3-10) 
P s-—»~o ntl 


Similar equations can be written for a type one and type two system. 
However, they will not be considered at this point due to the fact that 


they are of such a form that a linearization procedure to solve the system 


of four equations in four unknowns is not possible. On the other hand 
16 


the problem can still be solved using the iterative procedures of 


program FISDEP in Appendix D. 

The given set of specifications has been used to find equations 3-8, 

3-9 and 3-10. They constitute a set of three equations in four unknowns, 
therefore, another equation is required. This equation can be an 
auxiliary equation by means of which the designer can force a specific 
real root of the system to be placed at a determined location, allowing 
in this form that the given root specification will have at least a real 
part less than one of the real roots of the system (for a second order 
uncompensated transfer function this procedure will insure the dominance 
of the given root specification). Under this consideration the needed 
fourth equation is: 


n+l 


_ k 
F(s) = >, Gy (K,./K, .Z,P) Tx (4-17) 
k=o 


If the coefficients G, (Kf, Z, P, Kt) of equations 3-8, 3-9, 3-10 and 


3-11 are expanded, the following matrix equation can be written: 


jou Sie Oia diy ry] c | 
Foy So2 403 Seg ree 
= (B-12) 
Sai “a9 S39 “sa "3 73 
S41 S42 S43 Aga) 4 | “4 
where: 
r, = KK Z 
Tr, =P 
r, = KK 
r= KK + K KZ 
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and ait are the constant coefficients of r, in equation i (i = 3-8,.etc) 


c, are the constant terms in equation i (i = 3-8,.......3-11) 

Equation 3-12 is a linear algebraic matrix equation of the form D.r=C. | 
Solution of this equation is r = p . C. If the found solution is used 

to determine the values of the unknown parameters the following two 
alternatives result: 

Key = (1 /2K ) (t+ (c4°-4r,r4) 


Ky = %3/Key 
2 =1y/K Key 


Pp =T, 


alternative two « 
K,. = €4/2K) (1 + 2-41) 172) 
Keg = %5/Kep 
a=) /5 Boo 
P= 7) 
The proper choice of either alternative depends mainly on the numbers 
involved and secondly on the designer criterion when both solutions 
are physically possible. On the other hand it might be possible that 
neither solution one nor solution two is physically possible or that the 
system becomes unstable with one or the other solution. In both cases - 
another combination of compensator devices might be used. 
All the development has been done so far considering an external 


loop gain fixed; however, some problems might require trading a variable 


external gain, Ko by a fixed transfer function gain, Ke Obviously if 
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this is the case the above development still holds and the solution of 


the problem gives as before two alternatives, in which the compensator 
pole and the tachometer gain have the same numerical value while two 
values for external loop gain are found at the expense of the compensator 
zero. The following solutions can be written: 


alternative one 


Ko4 = ,/2K) + 7-42 5,) 77) 
bis tr /K Ky 
K, = r/K, 
ae 
“ alternative two 


Ky = ,/2K) (1 - 2-402) 


. ae deere 
Z,=1)/K Ko, 
K =r. 
1 = 83/ F 
P= To 
Sgnsidadad-C... ane, SS ec opment Ge tie 
nsiaering a an a2 g a. 7 pme a 


foregoing case (G_, (s) = 1) can be done and the problem can again be 


c2 


solved by finding a solution to equation 3-12 with the following definitions: 


Th) = KK A 
: ry =P 
ar: 
r= KOK, + KEK P 


Unlike the case where the filter is inside the internal loop, a unique 


solution is obtained in this case. On the other hand whether the external 
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gain or the transfer function gain is the variable, the equations con- 


taining the unknown parameters have the same form, that is 


P= ro 

K. = C, - TP)/K. 
K, = r./K, 

Z= r /KK, 


A computer program capable of handling all the cases considered was 
written by the author and is shown in Appendix D. 

As Examples 

Several examples will show the feasibility of using the 

computer program of Appendix D for different uncompensated transfer 
functions. 
Example 3.1 
Design a filter compensator and find values for forward and feedback 
gain which will give the system shown in Figure 3-2 a constant error 


coefficient of 2.0 and a pair of complex roots located at g = 0.5 and 


W_=5.0. It is also specified that a real root be located at = ==0 0. 


_ ki 
(S+5) (S+7) (S+10) (S+12) (S+15) 


Figure 3-2 


Block Diagram for Example 3.1 


20 


INPLT DATA 


ee ee a ee 


SPECIFICATION FOR ZETA [Sccavcecccescecvccccvccccve NeSONOGCE 95 
SPECIFICATICN FOR NATURAL FREQUENCY [Seccccesccces 9e-S50900E 91 
SPECIFICATION FOR REAL RCCT [Sccescceseccvcccceee ~“MeBCHHOE OL 
SPECIFICATICN FOR ERRCR CCEFFICIENT [Ssccesevcvees 0.e20CGC0E 31 
SPECIFICATION FOR EXTERNAL LOOP GAIN ISccccccncce GeS5N000EF O1 


COEFFICIENTS OF NUMERATOR POLYNOMIAL IN DESCENDING CRDER 
Qe-1CCCOE C1 

CGEFFICIENTS OF DENOMINATOR POLYNOMIAL IN DESCENDING ORDER 
DelCOOCE O1 0.49000E G2 Ne 9290CE 03 92 84950E 04 
N.37350E 05 9e63090CE 95 


REAL ZERCS GF UNCOMPENSATED TRANSFER FUNCTION 


920 


CCMPLEX CONJUGATES ZEROS OF UNCOMPENSATED TRANSFER FUNCTION 


9.0 0." 


REAL POLES OF THE UNCOGMPENSATED TRANSFER FUNCTION 
~O0-5CCOCE (1 ~O.79N0CE 1 ~9.10900E G2 ~3.120C0E 02 


~TelLECOCE 2 


CCOMFLEX CONJUGATES POLES OF UNCOMPENSATED TRANSFER FUNCTION 


9.0 0.¢ 


a1 


FILTER CUTSIDE INTERNAL LOOP, 


FORWARD GAINeccecccccceeeeK= 
TACHCMETER GAINeccecccces KT 
COMPENSATOR POLEccecccceePC= 
COPPENSATGR ZEROcccccceeelC= 


CCMPUTED ERROR COEFFICIENT IS 


COEFFICIENTS OF THE C.E. 


N0e27782E 


0.986C6E 


0.-10C00E 


06 


04 


Ol 


THE CCMPENSATED 


REAL PART 


-0.47188E 


-0-8COO1E 


-0.16376E 


-0.16376E 


~0.25000€ 


-9.25000E 


01 


Ol 


02 


02 


ol 


ol 


O-18007E 06 


O-1LOO1OE 04 


EXTERNAL GAIN SPECIFIED 


Oe1L0995E 05 
0244516E 90 
0-14699E O01 
0-33689E 91 


0-20009E Ol 


IN ASCENDING ORDER 


0e54732E 05 


9-50470E 02 


SYSTEM ROOTS ARE 


IPAGINARY 
0.0 
0.9 
~O-SLI9STE 01 
O-5LI197E O1 
-0.43301E O1 


0.43301E O1 
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Example 3-2 


2.0(S+5) (S+8) (S+10) (S+3~-5}) (S+3+5)) 
(S+3 . 5) (S+6) (S+9) (S+15) (S+1.5+3j) (S+1.5-3)) 


Figure 3-3 
Block Diagram for Example 3-2 


Design a filter compensator and find values for external loop and feedback gains 
which will give the system shown in Figure 3-3 a constant error coefficient of 5.0 anda 
pair of complex roots located at 7 = 0.5 and my = 5.0. Itis also specified that a real 


root be located at g, =-10.0. 


SPECIFICATICN FOR ZETA [Sec cc cccccccccscescencece GeS5ONCOE wt 


SPECIFICATION FOR NATURAL FREQUENCY [Scecccccccce 9e50000E 01 


SPECIFICATION FOR REAL ROOT 


TScccccccccccccccccecs ~O0.-100C0OE 02 


SPECIFICATION FOR ERROR COEFFICIENT ESecccccccsece O-50000E 91 


SPECIFICATION FOR TRANSFER GAIN ISeccoccececcecce 9.20000E 02 


COEFFICIENTS OF NUMERATOR POLYNOMIAL IN DESCENDING ORDER 


§:aisoce G4 8 £3888E 


COEFFICIENTS OF DENOMINATOR 
O-1CCOQOE O1 0. 36500E 


9-12514E O5 0.28603E 05 


REAL ZEROS OF UNCOMPENSATED 


-0.5C0C0E O1 -0.8000CE 


CCMPLEX CONJUGATES ZEROS OF 


-0-30000E 01 0. 5000CE 


E 92 9¢34200E 03 02.22020E 04 


POLYNOMIAL IN DESCENDING ORDER 
02 0-4957T5E 03 0233154E 04 
Oe31894E 05 


TRANSFER FUNCTION 


01 ~-O-1ND00E 02 


UNCOMPENS ATED TRANSFER FUNCTION 


01 


REAL POLES OF THE UNCOMPENSATED TRANSFER FUNCTION 


—-0.35CO00E O01 -0.60000E 
CCMPLEX CONJUGATES POLES OF 


~-O.15000E 01 £-0-30000E 


01 -9-90000E O01 -0.15000E 02 


UNCOMPENSATED TRANSFER FUNCTION 


O01 
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FILTER OUTSIDE INTERNAL LOOP, TRANSFER GAIN SPECIFIED 


EXTERNAL GAINewccccccces ce ck= %-7T4439E 40 
TACHOMETER GAINecocvccce cKT= Oe TET5STE d1 
CCMPENSATOR POLEceeccceeelFC= M.1L1L648E 90 
COMPENSATOR ZEROcccccceeel(= Oe 7TE8TS5S7TE IL 
CCMPUTED ERROR COEFFICIENT IS = Q.500C9E O1 


“COEFFICIENTS OF THE CoE. IN ASCENDING ORDER 
Oe1LSI3Z6E 67 02-63472E 06 Qe 1lS5374E 46 


0-45667E 05 Ne82727E 04 2e¢B86074E 03 
0-61387E 02 0233295E 01 


THE CCMFENSATED SYSTEM ROOTS ARE 

REAL PART IMAGINARY 

~Q.51883E O1 9.0 

QO-LLT6CGE O1 —-0-40260E 01 

O-LI76CE G1 0.40260E O1 

“O5e64411E O1 ~“O31541E O01 

~O.64411E C1 NOe31541E C1 

-0.213594E 01 “Q.10980E 02 


~OUel3594E O01 Ce10980E 02 
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Example 3-3 


K(S+90) 
(S+8) (S+15) (S+20) (S+25) (S+30) (S+50) (S+60) 


Figure 3-4 
Block Diagram for Example 3-3 


Design a filter compensator and find values for forward and feedback gains which 
will give the system shown in Figure 3-4 a constant error coefficient of 0.5 and a pair 
of complex roots located at $ = 0.5 and We = 20.0. Itis also specified that:a real 


root be located at en = 100i, 


INPLT DATA 


SPECIFICATION FOR ZETA [Sccccccccceccccccccecccce 0-509600E 09 
SPECIFICATION FOR NATURAL FREQUENCY [Scccccccccce 0-20000E 22 
SPECIE RCA TON FOR REAL RDOT [Sccccncccccccccccces ~CelOONDE 03 
SPECIFICATION FOR ERROR COEFFICIENT [Scccccccccne 0.-50009E 09 
SPECIFICATION FOR EXTERNAL LOOP GAIN TScccccccone 0.30000E 01 


CCEFFICIENTS OF NUMERATOR POLYNOMIAL IN DESCENDING ORDER 
Oe1OCNCE O1 0.-SOQ0CE 02 
COEFFICIENTS OF DENOMINATOR POLYNOMIAL IN DESCENDING ORDER 


OC.1CQCOE O1 0-20800E 03 NeLT475E 05 Oe 76700E 0 
0.18972E 08 0-.26382E 09 918990E 10 DQeS54Q0Q00E 1 


6 
6) 
REAL ZERGS OF UNCCOMPENSATED TRANSFER FUNCTION 


-0.SCOQ0E 02 


COMPLEX CONJUGATES ZEROS CF UNCOMPENSATED TRANSFER FUNCTION 


O.€ 0.0 


REAL FCLES OF THE UNCOMPENSATED TRANSFER FUNCTION 
-0-8CO0CCE O01 -0.15000E 02 -9-2000CE 02 “9e25NO0E 92 


-0.-3C00CE C2 -~0.50000E 02 ~2e60000E 92 


CCMPLEX CONJUGATES POLES CF UNCOMPENS ATED TRANSFER FUNCTION 


0.0 0.0 
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FILTER INSIDE INTERNAL LOOP, EXTERNAL GAIN SPECIFIED 


CHCICE ONE 


FORWARD GAINeccccccceccesekK= Gel2230E 99 
TACHOMETER GAINecccccceeeKlT= 0.26480E 90 
COMPENSATOR POLEcececvccese oFC= 0-99600E 02 
CCMPENSATOR ZERO ce cccccceell= 02e81438E 01 
COMPUTED ERROR COEFFICIENT IS = 0.50000E 00 


COEFFICIENTS OF THE CeEe IN ASCENDING ORDER 
0.80676E 12 0.25429E 12 Oe31IT21E 11 


0e21859E 10 0.95365E 98 0.25075E O7 


0.38192E 05 0-3076NE 03 O-10000E 01 
TRE CCMPENSATED SYSTEM ROOTS ARE 
REAL PART IMAGINARY 

-0.8C843E Ol 0.0 

~O-11758E 02 0.0 

~0-1CCO0E O02 0.17320E 02 

-Q.10000E 02 -O0.17320E 02 

-0.71132E 02 0-0 

~0.48312E 02 0-25473E 02 

-0.48312E 02 ~0.25473E 02 


-0.99999E 02 9.20 
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CHOICE TWO 


FCRWARC GAIiNeccccccccceceeK= COeB87914E 38 
TACHCMETER GAINecocccvceeKT= 0.-36838E 96 
CCMPEARSATOR POLEcccccvceePC= O.99609E 02 
CCMPENSATOR ZEROcccecccecell= O.11L329E 92 
CCMPUTED ERROR COEFFICIENT IS = 0-500G)E Of 


COEFFICIENTS OF THE CeEe IN ASCENDING ORDER 


0. 8C676E 


0+21859E 


0.38192E 


IZ 


10 


05 


THE CCPFENSATED 


REAL PART 


~}. 80843E 


-O.11758E 


~Q. 1GGQ0E 


~0.1CGOC0E 


-O0.71132E 


~-0.48312E 


-0.46312E 


-0.99999E 


eR | 


02 


02 


02 


02 


G2 


02 


C2 


0.25429E 12 Oe31721E 11 


0.95365E 08 9.25075E O7 


9.30760E 03 9-10000E 01 


SYSTEM ROOTS ARE 


IMAGINARY 


0.6 


0.0 


O0.17323E 02 


-9.17320E 02 


020 


0.25473E 02 


-0.25473E 02 


9-0 


od 


0€ 


Example 3-4 
af 


5000,000 (S+2+104) (S+2-10}) 


(S+12) (S+14) (S+8-8j) (S+8+8)j) (S+4 ..5-5j) (S+4.5+5j) 


Figure 3-5 
Block Diagram for Example 3-4 


Design a filter compensator and find values for external and feedback gain which 
will give the system shown in Figure 3-5 a constant error coefficient of 0.2 and a pair 
of complex roots located at § = 0.5 and We =4,0. It is also specified that a real root 


be located at c. = =15 0% 


INPUT CATA 


— me wo ee ee wee 


SPECIFICATION 
SPECIFICATION 
SPECIFICATION 
SPECIFICATION 
SPECIFICATION 


COEFFICIENTS 
0. 1CGQ0GE 
COEFFICIENTS 


0. 1C000E 
9. 48636E 


REAL ZEROS OF 


0.0 


FOR ZETA TSeccccecccn ce ccverssccesece % e500G0E 


FOR NATURAL FREQUENCY [Secccccccece 2e40000E 


FoR REAL ROOT TScccccccccccccccccce ~O.1L500G0E 


FOR ERROR COEFFICIENT [Seccccvcccce Je 2Q0Q0C0E 


FOR TRANSFER GAIN IScccocccsccccccce 0. 5C000E 


OF NUMERATOR POLYNOMIAL IN DESCENDING CRDER 


Ol 0.-4000CE O1 VeAALO4O0E 03 

OF DENOMINATCR PCLYNOMIAL IN DESCENDING ORDER 
Cl 0-69D00E 02 0e21522E 04 Q2.39970E 05 
C6 0.409094E 07 9e22356E 08 YeBLICGE 08 
09 O.15764E OS 


UNCOMPENSATED TRANSFER FUNCTION 


CCWPLEX CONJUGATES ZEROS CF UNCOMPENSATEO TRANSFER FUNCTION 


-Q.2CG600E 


REAL PCLES OF 


~4. 30000E 


~0.14Q0CE 


ol Q-1N000E 02 


THE UNCCMPENSATED TRANSFER FUNCTION 


01 -O.60000E O1 -929G000E 01 ~Je120GCE 02 


G2 


CCOMFLEX CONJUGATES POLES CF UNCOMPENSATED TRANSFER FUNCTION 


~0.8CG0C00E 


Ol GO.-8N9D0E 01 —9e45000E 91 %e 50000E 1 
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eae) 


Ol 


2 


on 


6 


FILTER INSIDE 


CHCICE CNE 


EXTERNAL GAINecoccccccceeckK= 
TACHOMETER GAINeccocccceeeKT= 
CCMPERSATOR POLEccccccceePC= 
COMPENSATOR ZEROcccccceeelZl= 


CCMPUTEC ERROR COEFFICIENT IS = 


0-18808E 91 
Oe1LL514E 00 
0.-58995E 91 
O-19018E OL 


0-20000E 00 


COEFFICIENTS OF THE CeE~ IN ASCENDING ORDER 


OellL16CE 10 
0-21427E 09 
0.72216E 06 


0-74899E 02 


REAL PART 
-0-22132E O1 
-0.13493E 02 
-0.20000E O1 
“0247320E 01 


~0.86161E 01 


0.12874E 10 


0.46067E 08 


0-52667E 05 


QO-1O0000E 01 


THE CCMPENSATED SYSTEM ROOTS ARE 


IMAGINARY 


0.0 


-0.38669E 01 


"0 .-34641E 01 


0 .69353E 01 


-0.87766E Ol 
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9e66220E 09 


90-68786E 07 


90225593E 04 


REAL PART 


-0-15003E 02 


-0213493E 02 


—e20000E 01 


-9.-47320E 01 


“Oe86161E O1 


INTERNAL LOOP, TRANSFER GAIN SPECIFIED 


IMAGINARY 


3-0 


0238669E Ol 


0«34641E OL 


0.69353E O01 


Oe877T66E O1 


CHOICE TWO 


EXTERNAL GAINeceocccccccvceekK= 
TACHCWETER GAINecoccccce ckKT= 
COMPENSATOR POLEcces cccce cofC= 


COEFFICIENTS OF 


O-11160E 1¢ 


0e21427E 0S 


C.7Z2216E G6 


C-74899E 02 


THE COMPENSATED 


REAL PART 


-0e22132E O1 


-0213493E 02 


-~0.2C00CE O1 


~“0047320E O1 


-0.-86161E O1 


COMPENSATOR ZEROccvccccesll= 


COMPUTED ERROR COEFFICIENT IS 


THE CLE. 


90-12874E 10 


0-46067E 98 


Oe526E67E OF 


O-10000E 01 


SYSTEM ROCTS ARE 


IMAGINARY 


0.0 


~0-38669E O1 


~0.34641E O1 


~0 -69353E Ol 


—O0.87766E O1 
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O.218S7E 00 
O.11514E 00 
0.58995E O1 
0-16335E 02 


0. 20000E 00 


IN ASCENDING ORDER 


%e66220E 09 


Q-68786E D7 


9e25593E 04 


REAL PART 


-0.15003E 02 


~9213493E 02 


-9.20GO00E 01 


-0.47320E O1 


“OQ -B6161E OL 


IMAGINARY 


0.9 


0.38669E O01 


0.34641E 01 


9.69353E Ol 


0.87766E Ol 


B. ROOTS AND SENSITIVITY SPECIFICATIONS : 


The problem of design, as was stated in section A, requires as 
final step the selection by the designer of the "best" system if more than 
one satisfies the specifications. The "best" system cannot be defined 
by a simple statement. A given set of specifications can give one or 
more solutions, and if all of them are physically possible the selection 
of the "best" system probably will be made taking into account cost, 
availability of materials, etc. 

It is possible for the designer to choose a system in which the 
dominant root, for example, satisfies the required time response but 
small variations of a system parameter change its location to positions 
in which it is not dominant anymore or makes the system unstable. 
Therefore, it seems obvious that a knowledge of these variations will 
help the designer choose the "best" system. 

One of the many definitions of sensitivity is [4]: The sensitivity 
of P with respect to Q is the percentage change in P divided by that 
percentage change in Q which causes the change in P, with all changes 


considered differentially small. 


Q d(InQ) dQ/Q 
Kokotovic and Siljak [5] developed algebraic equations to apply i 


Ae a(QinP) _ GE/P _ 2 (aP/aQ) (3-13) 


the concept of sensitivity as a measure of the quality of the system. 
They used the microscopic definition of equation 3-13 and the general- 
ized Mitrovic equations, equations 2-3 and 2-4, to develop algebraic 


. sensitivity equations of a system dominant roots with respect to a 
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SSS SS Se ae ne 


system parameter. If the parameter considered is the system gain, Kf, 


the equations can be written as follows: 


S F = = Re (3-14) 
Kf €D 
Ww 
See ~K, Dee (3-15) 
D 
where 
a = a Fj C B 
a D — pwn = it 1 
Kf Kf 
A, Bo A, My Co Pe 
n n 
7 k k _ - k 
A= (-1) "KAW “UL _4( ¥) A, i (-1) "KAW ULC ¥) 
k=o k=o 
n n 
_ k Kens 7 er : aan 
B= GD AW Uy ae) a (IAW UL CF) 
k=o k=o 
a= 5 ey DA 
, — Ew ky (4) a k OA 
k=0 OK, n k-l oF, = (-1) OK, Ww ULC F ) 
k=o 


and UL ( $), UL ( ¥) are the derivatives of UL Ff) and U 


are related by the equation 


UGS ea GU a) Ud Oh) 
Equations 3-14 and 3-15 constitute the required equations to be used in 
a design problem in which root sensitivity is a specification. They have 
to be used in conjunction with root specification, that is, root sensitivity 


specification is not independent of roots specification. 
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i. General Statement of the Problem 

Design the adequate compensator devices required to give 
(for a specific transfer function) a closed loop time response determined 
by root locations at ¥ = ¢, . i we = We It is further required that changes 
in system gain by X percentage should not change ¥ by more than Y 
percentage. 

The statement of the problem implies three specific ations, 
namely damping ratio, natural frequency and sensitivity. Therefore, 
considering the system gain as a parameter, compensator devices with 
two unknown parameters must be inserted into the system in order to 
satisfy the specifications. Then the application of equations 2-3, 2-4 
and 3-14 give a set of three nonlinear algebraic equations in three | 
unknowns, as required for solution. Several combinations of compen- . 
sator devices in cascade and/or feedback may be used. Typical 
combinations are acceleration and tachometer feedback, or one or more 
identical sections of filter in cascade. If acceleration and tachometer 
feedback are used the characteristic equation of an n order uncompen- 
sated system (assuming n =2) has the following form: 


n 
? 2 —- - 


k=o 


Replacing 3-16 in 2-3, 2-4 and 3-14 with 


- - 2h < = 
See XY Constant S. 
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the following equations can be written 


K.D 
=. ft KE 7 

ae =D G17) 
. k k 
s (-1)"G, (K,, KK JWI UL _1( $.) = 0 (3-18) 
k=o 
= k k 
>. (-1) "GK, KK WU, ( eS )j= 9 (3-19) 
k=o 


Analytical solution of equations 3-17, 3-18 and 3-19 by hand is, if not 
impossible (it can be done for very simple systems) at least difficult 
and time consuming, but if a digital computer is available the numerical 
solution is fast and accurate. 

A digital computer program capable of handling this problem 
is presented by the author in Appendix D. 
Example 3.5 
Determine the forward, acceleration and tachometer gains of the system 
shown in Figure 3-6 such that the closed loop system has a pair of poles 
at <, = 0.5 and Ws = 5.0. Itis also required that if variations of gain 


are 10% the corresponding of o being 5%. 


Kf 


R(s)+.. + + C(s) 
= is S(S+1.5) (S+4.5) | 
| aE ae eee en K ¢ Se 
j 
i 
| 
sie: ——{ KS }—— | 
eee. eee we 
Figure 3.6 


Block Diagram for Examples 3.5 and 3.6 
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In order to be able to use the computer program some preliminary calcu- 
lations are required. 


The open loop transfer function can be written 


G(s) S 
G_(s) = a ae i = of as dt a oo 
a 1+G(s) K $°+KS S°+(6+K K )S°+(6.75+K.K)S 
a ie fa £t 
The closed loop transfer function is 
G_(s) K 
f 
C(s) 2 = = 
R(s) 1+G (s) 3+ (6+K.K_)S7+(6.75+K K.)S+K 
fe) fa ra f 


Therefore the characteristic equation can be written as 


CE(s) = A(4)S° + A(3)S7 + A(2)S + A(1) 
where 
A(4) = 1.0 


A(3) = G.0+KK =6.0+XZ 


f 
= 6.75 +KK, = 6.75 +X 
A(2) = 6.75 + KK, = 6.75 + XY 


A(i) = <7 * 


taking the partial derivatives of the above coefficients w.r.t.X (K,) 


AD(4) = 0.0 
AD(3) =Z 
AD (2) =Y 
AD(1) = 1.0 


These are the basic calculations required as input to the program in 
Appendix D, however, the details concerning its use are explained in 


the program itself. 
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Two solutions to this problem are shown in the following pages. 


In both solutions the specifications are exactly satisfied. However, 

in the first one the real root that the designer does not want to be 
dominant is located at 7, = -4.50 while in the second one that root 

is located at C.. = -21.25. Therefore, it is obvious that between these 
solutions (the system may have other solutions) solution #2 should be 
better. It is worthwhile to note that the initial values for the unknown 
Parameters are quite different from the actual ones. However, a good 


guess implies less iterations and therefore less time of calculation. 


ag 


INITIAL VALUES FOR UNKNOWN PARAMETERS ARE 


0.69000D C3 
0.30000D CC 


SPECIFICATICN FOR 
SPECIFICATION FOR 
SPECIFICATION FCR 
SPECIFICATION FOR 


SPECIFICATION FOR 


0.240000 00 0-400000-01 


BANDWIOTH [5S 


ERROR COEFFICIENT IS 


ROOT SENSITIVITY IS 


DAMPING RATIO IS 


NATURAL FREQUENCY IS 


ITeUNKNOWN PARAMETERS 


FORWARD GAIN IS 


TACHOMETER GAIN [5S 


ACCELERATION GAIN 


NONE 


NONE 


0.20000D 00 


0.-50000D 90 


0.500090D 01 


0.531250 063 


0-23435D 00 


Is 0.238118D-01 


NUMBER OF ITERATIONS IS 


40 


65 


TI.CALCULATED SPECIFICATIONS 


SENSITIVITY IS 0-20000D 900 
ERROR COEFFICIENT IS 0640476D 01 


ITTeSTABILITY INFORMATION 


COEFFICIENTS GF THE CoE. IN ASCENDING ORDER 


0.53125D C3 0.13125D 03 0.26250D 02 


C.-10000D C1 


THE COMPENSATED SYSTEM ROOTS ARE 
REAL PART IMAGINARY 
~C.25000D Cl —-0243301D 91 
-0.250000 Cl Q.43301D 01 


—-0.21250D 02 0.0 
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INITIAL VALUES FOR UNKNOWN PARAMETERS ARE 
0-10000D G3 0.300000 Gd 0.300000-01 
C.20000D 60 


SPECIFICATICN FOR BANDWIDTH IS NONE 
SPECIFICATION FCR ERROR COEFFICIENT IS NONE 
SPECIFICATION FCR ROOT SENSITIVITY IS 0.200000 90 
SPECIFICATION FOR DAMPING RATIO IS 0.50C00D 00 


SPECIFICATION FOR NATURAL FREQUENCY IS 0.500000' 01 


TeUNKNOWN PARAMETERS 


FORWARD GAIN I$ 0-11250D0 03 
TACHOMETER GAIN {S$ 02362220 00 
ACCELERATION GAIN IS 0.31111D-01 
NUMBER OF ITERATIONS IS 13 
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II CALCULA TED SPECIFICATIONS 


SENSITIVITY IS 0.200000 00 


ERROR COEFFICIENT IS 0.236840 01 


ILTT~STABILITY INFORMATION 


COEFFICYENTS CF THE C.F. IN ASCENDING ORDER 


02112500 C3 0.475000 02 0.950000 01 


€.100C0D C1 


THE COMPENSATED SYSTEM ROOTS ARE 


REAL PART IMAGINARY 


-G.45000D C1 0.0 


-0.2500CD Cl ~0.43301D O01 


-€.250000 01 0.433010 O01 
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Example 3.6 


Redesign the system of Figure 3-6. Use in this case a root sensitivity 
constant such that variations in gain of 10% produces variation in ¢ 
of only 1% (Se = 0.1). 

As in example 3.5, two solutions are presented. Comparing these 
solutions with those of example 3.5 it should be noted that as the system 
gain is increased the sensitivity constant is reduced and the required 
roots are more dominant (solution two of example 3.6). This was also 
verified working the problem with sensitivity specification Se = 0.3 


in whose caSe the resultant parameters for two solutions were: 


Solution One Solution Two 

K = 133.59 K = 337.24 

K, = 0.3366 K. = 0.2540 ’ 
K =0.325 K = 0.0370 

a a 

Oy, = 75-343 G,, = 713.49 
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INT TIAL VALUES FCR UNKNCWN PARAMETERS ARE 
C.8090GC 03 Ge 2CO9UD O01) Je 30000 0-901 
Ceo30utlel Cf 


SPECIFICATICN FCR BANDWICTH IS NONE 
SPECIFICATICA FCR ERRCR COEFFICIENT IS NONE 
SPECIFICATICA FCR ROOT SENSITIVITY IS VelLIWOGID gs 
SPECIFICATICA FCR CAMPING RATIO IS 9e50G0G970 9) 
SPECIFICATICAN FCR AATURAL FREQUENCY 1s 5-50009D a1 


FCRWARD GAIA I§ J.97T57ID 92 
TACHOMETER GAIN IS Ge38704D 06 
ACCELERATICA GAIA IS 92297510-G1 
NUMBER OCF ITERATICAS IS ae 
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TTCALCULATED SPECIFICATICNS 


SENSITIVITY IS 92109060 00 
ERRCR CCEFFICIENT IS Is21919D O1 


TTTeSTABILITY INFCRMATICN 


COEFFICIENTS CF THE C.sE. IN ASCENDING ORDER 


CeS7E7T10 O02 02445140 62 3289028D O01 


Ceo1NGICe Cl 


THE COMPENSATED SYSTEM ROOTS ARE 
REAL PART IMAGINARY 
~U2e39028D 01 0.6 
-0.250GCD O1 7024336010 O1 


-0.2500CD 01 0.43301D CL 
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INITIAL VALUES FCR UNKNGWN PARAMETERS ARE 


ce6000CCD €23 


£.3000CC 


SPECTIFICATICN 


SPECIFICATICN 


SPECIFICATICNA 


SPECIFICATICN 


SPECIFICATICA 


cc 


FCR 


FCR 


FCR 


FOR 


FCR 


026246500 uv) 0e49C0UD-01 


BANDWICTH IS 


ERRCR CCEFFICIENT IS 


ROOT SENSITIVITY IS 


CAMPING RATIC IS 


NATURAL FREQUENCY IS 


Sseaeweenes f2oew ee eee 


FCRWARD GAIA I§ 


TACHOMETER GAIN IS 


ACCELERATICA GAIN IS 


NONE 


NONE 


Ge1l0C00D Be 


O.FOCG90 00 


02500000 21 


0210649D 04 


0e21714D 09 


NUMBER OF ITERATICAS IS 


47 


02.39061D-01 


25 


ITe-CALCULATEC SFECIFICATIONS 


SENSITIVITY IS 


ERRCR COEFFICIENT IS 


UelOVOGD OO 


0.44748D 01 


ITIL. STABILITY INFCRMATION 


ee ee a a ee 


COEFFICIEATS OF 


C.10649D 04 


f.1006CD O1 


THE COMPENSATEC 


REAL FART 


~U.25000D C1 


-9.250CCC Cl 


-0.42597D €2 


on Tee eae ee eee 


THE CeEe IN ASCENDING CRDER 


02237990 u3 0.475970 02 


SYSTEM ROCTS ARE 


IMAGINARY 


0.433610 61 


0.433010 Ol 


0.0 
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IV. FOUR SPECIFICATIONS DESIGN 


The frequency response specification which is most often specified 


is the bandwidth. It can be defined as the angular frequency, W,, at 


b 
which the magnitude of the output to input (closed loop transfer function) 
has suffered an attenuation M (often expressed in decibels). If the 

attenuation is three decibels (.707) then the bandwidth is usually called 


the "3 db. bandwidth." (In this paper this will be the bandwidth to be 


used.) 


A. ROOTS SENSITIVITY AND BANDWIDTH SPECIFICATIONS 

In Section III the significance of the root specifications as a 
measure of the time response of the system and the sensitivity as a 
measure of the quality of the system were stated. If bandwidth is now 
added as a specification then it can be said [4] that it is significant 
because it indicates rise time or speed of response, it measures in 
part the ability of the system to reproduce the input signal, and it 
approximately describes the filtering characteristics of the system. 

l, Statement of the Problem 

Design appropriate compensator devices to give (fora 

specific transfer function) root locations at ae = , vi = Was 
It is also required that the compensated system have a bandwidth We =B 
and that variations in gain of X percentage produce variation in ., 


of only Y percentage. 


AQ 


The problem as stated here is the same as Section III with 
the additional constraint of bandwidth. That is, in addition to equations 
2-3, 2-4 and 3-17, a fourth equation must be found. 

From the definition of 3 db bandwidth the following relation- 


ship can be written: 


Cc (jw, ) 
b 1/2 
ee Cee = = (1/2 zi 
RGW ) M = (1/2) (4-1) 
b 

If the system is to have a bandwidth Wy = B then equation 4-1 can be 

written as 
CGB) = 0.707 (4-2) 


R (jB) 

If both sides of equation 4-2 are squared then it can be expressed as a 
polynomial in B whose coefficients are functions of the unknown 
parameters, that is, if a closed loop transfer function is assumed to 


be of the form 


m m-l 
+ eoeeeoes eee + 
ef. Se eee ot) 
R(s) A,s" Pe aia ee +A +1 
then 
10) ..,m-l 
CUB) _ oo + Q, GB) a ae 
R(jB) - Le 
A, (iB) + A, (jB) Heese tA 
assuming m and n to be even 
-2 2 m-l 2 
Gp} 2 +Q.B" ee 
CB q (1Q,B £Q,8 Fee tQ yt ¢Q,B Hi. ..+Q_B) 7 
R(jB - “ oon- 7 
GB) (HA, BTA, BTS. 4A » + (A,B Ce a 
nh 


n+ 
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or 


gnol 2 


0 = GAB +...4A +...+A 8) 


— I n+] (4-3) 


2 m-l 


m 2 
= 2l+O,8 #.+.Q 47) a an +...+Q_B) ] 


Equation 4-3 is a rather complicated nonlinear algebraic equation which 
must be combined with equations 2-3, 2-4 and 3-17 to give the required 
set of four equations. Compensator devices must now be added to the 
system in order to satisfy the specifications, that is they altogether 

must introduce three parameters to combine with the system gain and 

form a four parameters design problem. 

Example 4,1 

Design a compensator filter and find the values for forward and tachometer 
gains to give the system shown in Figure 4-1 a 3 db bandwidth at 

W. = 6.0 radians/second, a pair of complex poles at $ = 0.5, Wo = 5.0 


b 


and a root sensitivity constant of 0.2 
Kf C(s) 


| S(S+1.5) (S+4.5) 


ee - 


Figure 4-1 
Block Diagram for Example 4.1 
Solution of this problem was obtained using the computer program of 


Appendix D. As in Section III the characteristic equation must be found: 
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The open loop transfer function is 


K, (S42) 

G(s) = G(s)G_(S) = gigyp) G1 .5) (844.5) + (S4P)K «5 

where 
Be 

G(s) = Sigri.5) (614.5) + KK, 

and 
oo) D2 

S. (s) = oo & 

the closed loop transfer function is 
K,(S+2) 
a a a ee 
R(s) S(S+P) (St1.5)(S+4.5) + (S+P)K KS + K, (S+Z) 


and the characteristic equation can be expressed as 


CE({s) = A(s)S" dy A(a)s® + A(3)8” + A(2)S + A(1) 


where 
A(1) = KZ 
A(2) = 6.75P + K KP Hf: Ky 
A(3) = 6.75 + 6P + KK. 
A(4) =6+P 
A(5) = 1.0 


Since root sensitivity has been specified the partial derivatives with 


respect to the gain K, are required 


f 


me 


With the variables A(i) and AD(i) G = 1,....5) and the corresponding 


specifications as input to the program, a solution is shown in the 
following pages. This is not a unique solution, it is just one of the 
many solutions the problem may have. Also, it is worthwhile to observe 
that the solution doesn't guarantee dominance of the given roots but it 
does guarantee that a pair of complex roots will be at the specified 
location and that the other specifications are also met. On the other 
hand, the fact that information of all the poles of the characteristic 
equation is obtained makes it possible to investigate by inspection 


whether or not the given roots are dominant. 
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INITIAL VALUES FOR UNKNCWN PARAMETERS ARE 
0.100000 Q3 0. 70000D 00 


0.300000 01 0.30000D 00 


SPECIFICATICN FCR ERROR COEFFICIENT IS 


SPECIFICATION FOR BANDWIDTH IS 


SPECIFICATION FCR ROOT SENSITIVITY IS 


SPECIFICATION FOR DAMPING RATIO IS 


SPECIFICATION FOR NATURAL FREQUENCY IS 


TeUNKNOWN PARAMETERS 


— ee ee es 


FORWARD GAIN I§ 


TACHOMETER GAIN IS 


COMPENSATOR PCLE IS 


COMPENSATOR ZERC IS 


NUMBER OF ITERATIONS I$ 
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0249027D 02 


0.916470 00 


0.389180 OL 


0.13040D 02 


34 


0.10000D 02 


NONE 


0.600000 O1 


0.200000 00 


9.500000 00 


0.500000 O01 


TIT.CALCULATED SPECIFICATIONS 


a w2e@2eaeese2esee2en2 = = 


ABS.VALUE OF CLTF AT B IS 0.707110 00 


SENSITIVITY IS 


0.200000 00 


ERROR COEFFICIENT IS 0.317860 01 


TITeSTABILITY INFORMATION 


COEFFICIENTS CF 


02639340 C3 


C.98918D O1 


THE COMPENSATED 


REAL PART 


-0.25000D 061 


-0.25000D 01 


-0.24459D0 Cl 


0.244590 G1 


THE C.E. IN ASCENOING CRDER 


0.250160 03 0.750330 02 


0.10000D 01 


SYSTEM ROOTS ARE 


IMAGINARY 


-0.433010 O01 


0.43301D 01 


02442620 G1 


0.442620 01 
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B. ROOTS, ERROR COEFFICIENT AND BANDWIDTH SPECIFICATIONS 


In earlier sections the basic concepts of bandwidth, error co- 
efficient, sensitivity and root. specifications have been discussed. 
However, since error coefficient was defined in Section III for a 
particular case the definition can now be generalized and the following 


equations can be written: 


K_ = lim G, (s) (4-4a) 
$70 
K ='lim SG (s) (4-4b) 
Vv ° 
$70 
; 2 
K =limS G (s) (4-4c) 
a ° 
$20 


where G(s) is the open loop transfer function of a negative unity feed- 
back system and the subscripts p, v and a refer to type zero, one or two 
system. 

The error coefficient can be interpreted in terms of the character- 
istics of system performance in both the time and the frequency domain [4]. 
From equations 4-4 it is observed that only one of the error constants 
possesses a finite, non-zero value for any given system and this value 
describes the low frequency behavior of the open loop transfer function 
(at low frequencies large error coefficients imply high gains). On the 
other hand the error coefficients can be interpreted in terms of the steady 
state errors; i.e., if E(s) is defined as the difference between the Laplace 
transform of the reference input signal (R(s)) and that of the output signal 
C(s) in a negative unit feedback system, and if E(s) possesses poles in 


the left half plane only then the steady state error can be written as 
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SR (s) (4-5) 
ss 1+G (s) 
$70 re) 


E = lim 


From equations 4-4 and 4-5 the following equations can be derived: 


1 _ 1, 
Eos =r (for R(s) = 3 ) (4-6a) 
p 
po= ob (for R(s)= 
cs k or R(s) = 52) (4-6b) 
Vv 

_ 1 _ tL 

Eos = x (for R(s) = 33) (4-6c) 


Therefore, K_, KY and K, are measures of the steady state errors if the 
input is a unit step, a unit ramp and a unit parabolic function 
respectively. 
Example 4.2 

° Design a filter compensator and find values for forward and feedback 
gains required to give the system shown in Figure 4-2 a pair of complex 
poles located at $ = 0.5,W, = Ww, = 9%. Itis also required that the 


system have a bandwidth (3 db bandwidth) at the frequency wy =6 


and a constant error coefficient K = 3.2. 


Figure 4-2 


Block Diagram of Example 4.2 
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The characteristic equation of the compensated system is 


E(s) = A(S)S* + A(4)S° +A(3)S" + A(2)S + A(1) 


where 
A(1) = 
A(2) =6.75P +K 
A(3) = 6P + 6.75 + KK IP 


A(4) =6.0+ KK P 


Likewise the bandwidth and error coefficient equations are 


i2 1/2 
[(K,Z) + (6 B)1} 
j<uay - gnatl’?-< ra 


[(a(s)B?-a(3)B°+A(1)]7 > tae A(4)B”) 3251/2 


K =K.Z/(A(2) - K,) 
Vv if 


t 
The acceleration feedback is probably not the most adequate device to 
be used in order to meet the specifications. However, the idea of the 
author was to have a case in which the nonlinearities introduced to the 
system of equations make difficult the convergence on the method to 
solve them. This is clearly shown in the solution on the next pages in 
which the specifications are met but the solution was obtained using 

a rather large number of iterations (661), see Appendix C, even though 
the termination criterion was chosen large enough (0.4). On the other 
hand the results can be considered good from the engineering point of 


view, 
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INITIAL VALUES FOR UNKNOWN PARAMETERS ARE 


C€.40000D ¢3 
0.-100C0D C2 


SPECIFICATION FOR 


SPECIFICATICN FCR 


SPECIFICATICN FOR 


SPECIFICATION FOR 


SPECIFICATION FOR 


FORWARD GAIN IS 


0.50000D-02 0.700000 0C 


0.40000D 00 


SENSITIVITY IS 


ERRCR COEFFICIENT IS 


BANDWIDTH IS 


DAMPING RATIO IS 


NATURAL FREQUENCY IS 


0.38524D 03 


ACCELERATION GAIN IS 0.69901D-02 
COMPENSATCR POLE I§ 0.129780 62 
COMPENSATOR ZERC IS 0.727680 00 
NUMBER OF ITERATICNS IS 661 
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NONE 


0.320000 01 


0.606000 91 


0.500000 09 


6.50C000 O1 


IIeCALCULATED SPECIFICATIONS 


ABS.VALUE OF CLTF AT 8 IS 0.70711D 00 
SENSITIVITY IS -0.35906D 00 
ERROR COEFFICIENT IS 0.32000D O01 


ITIeSTABILITY INFORMATION 


COEFFICIENTS OF THE C.E~ IN ASCENDING ORDER 


0.280330 063 0.47285D 03 0.119570 03 


0.216710 C2 0.100000 Oi 


THE COMPENSATED SYSTEM ROOTS ARE 
REAL PART IMAGINARY 
-0.70219D OC 0.0 
-0.250000 Ol -0.43301D Ol 
-0.250C0D Cl 0.43301D Ol 


-0.15969D 02 0.0 
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V. FIVE SPECIFICATIONS DESIGN 


In Sections III and IV the problem of design was solved assuming that 


a maximum of four specifications were required to meet. The problem 
however can be made much more complex if all the specifications 
studied so far are to be satisfied at the same time (this is not an un- 
realistic case but a desired one). 

It is easy to see that the system of equations which results when 
equations 2~3, 2-4, 3-17, 4-3 and 4-4 are combined is a very compli- 
cated one even though the system transfer function to be studied is 
the simplest one. This is due to the fact that since five specifications 
must be met the number of compensator devices has to be increased 
in such a way that five unknown parameters are introduced to the system. 
Under the assumption of variable transfer gain a possible block diagran 
of a system which is expected to meet five specifications is shown in 
Figure 5-1 (many combinations can be made but all of them must have 


a maximum of five unknowns including the system gain). 


Figure 5-1 


Block Diagram for Five Specifications Design 
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As in the examples studied in the foregoing sections this case can also 


be handled by the computer program included in Appendix D. The calcu- 
lations required to use the program can be summarized as follows: 
define 
2 
H(s)=K S +K5 
a t 


G(s 
1 + G(s)H(s) 


G(s) = 
then the open loop transfer function can be written 
G (s) = G_(s)G_(s) (5-1) 
fe) e c 


and the closed loop transfer function as 


C(s) G,{s) (5-2) 


R(s) 1+G (s) 
° 
From equation 5-1 the constant error coefficient is 


K = lim sic (s) (5-3) 
s70 


where 1 = 0, 1 or 2 corresponds to x = p, vor a for the three more 

common types of error coefficients. 

From equation 5-2 the characteristic equation can be written as 
CE(s) = 1+ G(s) (5-4) 

and the 3 db bandwidth equation as 

Cc Gw,) 


$$$ 


R (jW,) 


|G Ww ) 

s.. & 172 

2 = «S((1/2) (5-5) 
h + G Gw,)| 


Equations 5-3, 5-4 and 5-5 are the necessary equations required as 


input to the computer program. 
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The fact that five specifications must be met doesn't imply that 


the system can be forced by using a prescribed set of compensator 
devices to meet those specifications. Hence the designer must decide 
if all the specifications cannot be met, which of them should be 
sacrificed in favor of the others or in favor of the impossibility of using 
another type of compensator device. With this in mind and taking the 
block diagram of Figure 5-1 as an example, the flowchart of Figure 5-2 
was realized in the computer program of Appendix D. However, it is 
worthwhile to note that the program was written to handle the character- 
istic equation of the compensated system of three, four or five parameters 
and, therefore, the flowchart of Figure 5~2 is just an example that uses 
the configuration of Figure 5-1. 

Example 5.1 

Design a filter compensator and find values for forward and feedback 
gains required to give the system shown in Figure 5-3 a pair of complex 
poles located at & =0.5, ne = 5.0, a 3 db bandwidth at frequency 

B = 6.0 radians/second, a constant error coefficient KL =5.0 anda 


root sensitivity constant Se = 0.2. 


Kf 
S(S+1.5) ($+4.5) 


yremarnmentnnintnnneeey tenn Gayest NY eRe deme Ret int ia oni tm Weer aete ie aden 


Figure 5-3 


Block Diagram for Example 5.1 
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Read Specs. 
Wn, Wb,Kx, etc. 


Guess values of 
unknown params. 


Use subroutine 
to find coef. C.E 


fest 
Realizability peat 


and 
Stabilit 


Test 
Realizability 
and 


Test 
1 Realizability 
and 
Stability 


Figure 5-2 


Flowchart showing the design procedure (given a set 
of five specifications) for block diagram of Figure 5-1 
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Following the steps on page 62 equations 5-3, 5-4 and 5-5 can be 


written, for this particular case, as follows: 


K = K,2)/(A(2) - K,) 


f 
CUB) = (/y)!/2____Lana)-an@ye")?+aneyp-an ay) p72 
R(jB)} : 2 4 6 3 5.212 
[(A(1) -A(3)B° +A(5)B -A(7)B-) + (A(2)B-A(4)B° +A(6)B)“] 


CE(s) = a(5)S° + A(4)S? + A(3)S* + A(2)S + ACL) 
where 
A) = Kz 
=6.75P + K P+ 
A(2) 6.75 K _ Ke 
= + + 
A(3) 6P + KARE KK. 6.75 


A(4) =KK P+6 
fa 


the partial derivatives of the coefficients A(i) with respect to the 
system gain are given by 
AD(1) =Z 


AD (2) = K.P + 120 


AD(3) = K, 
AD(4) = 0.0 
AD(5) = 0.0 


With this data a solution was obtained and is shown in the following 
pages. Comparing this solution with examples 3.5 and 4.1, the following 
remarks can be written: 


a. In all the cases the specifications are exactly met. 
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INITIAL VALUES FOR UNKNOWN PARAMETERS ARE 


C-50000D 02 
C.390C0D C1 


SPECIFICATION FOR 


SPECIFICATION FOR 


SPECIFICATION FOR 


SPECIFICATION FOR 


SPECIFICATION FOR 


TeUNKNOWN PARAME TE 


—o§awe ae a ee 


FORWARD GAIN IS 


TACHOMETER GAIN IS 


G.90000D 00 9-13000D 02 
0.20000D-01 0.200000 00 


ROOT SENSITIVITY IS 
BANOWIDTH IS 

ERROR COEFFICIENT IS 
DAMPING RATIO IS 


NATURAL FREQUENCY IS 


RS 


0.511320 02 


0.87909D 00 


ACCELERATION GAIN IS 02120200-02 
COMPENSATOR PCLE I§ 0.38503D 01 
COMPENSATOR ZERO I§ ; 0212458D 02 
NUMBER OF ITERATIONS IS il 
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0.20G00D 00 


0.60009D O1 


0.320000 O1 


0.500000 900 


0.50000D 01 


IT CALCULATED SPECIFICATIONS 


ASBS.VALUE OF CLTIF AT BIS 0.767110 00 


SENSITIVITY IS 9-20N01D 00 


ERROR COEFFICIENT IS 0.32000D 61 


TITYeSTABILITY INFORMATION 


eee EET OS ee em ae ae om Re 


COEFFICIENTS CF THE CeEe IN ASCENDING ORDER 


C.636S9D 03 0.250190 03 0.75038D 02 


C.991L17D Cl C.10000D Ol 


THE COMPENSATED SYSTEM ROOTS ARE 


REAL PART IMAGINARY 


-0.250CCD O1 -0.43301D O1 


-€.25000D Cl 0.433010 91 


-0.24559D Cl -0.44100D C1 


-0.24559D Cl 0.441060 O61 
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Loy Dominance of the given pair of complex poles is only 
met, for the same system, with the three specifications problem. 

Ce System parameters for examples 4.1 and 5.1 are 
comparable. 

d. More than one solution can be obtained in each case 
(for simplicity in this paper only one solution is shown for examples 
4.l and 5.1). 

e. Each of the possible solutions is greatly dependent on 
the initial guess. 
Example 5.2 
The system shown in Figure 5-4 is completely unstable due to the 
presence of a pair of complex poles on the right hand plane. Design 
two identical sections of filter and find the forward and feedback gains 
required to give the compensated system a pair of complex poles located 
at F= 0.5, W_ = 5.0; a3 db bandwidth at a frequency W, = 1.0 


b 


radians/second; a constant error coefficient Ks = 2.0 and a root 


sensitivity constant Se = 0.2 


K(S+1), 
(S+2) (S+5) (S-2-5j) (S-2+5}) 


Figure 5-4 


Block Diagram for Example 5.2 and 5.3 
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Following the same steps as example 5.1 the equations can be written: 


K = (K.2")/290P" 
je) f 
CE(s) = a(7)s° + A(6)s° + A(5)S" + A(4)s° + A(3)87 + A(2)S + A(l) 
[(K,2)?+(K, p)241/2 
1/2 
= (1/2) = OO eee 
ROB) iB) [(A(S)B4-A(3)B°+A(1))~ + (A(2)B-A(4)B°) 2] 2” 2 


where the following definitions were made: 


AN(1) = ea 


AN (2) = 2KZ + KZ 


AN (3) = 2K 2 + K, 


AN (4) = 1.0 


A(1) = 290P* + KZ 


2 
A(2) = 580P + 163P + K KZ ae 2K 2 * Kz" 


2 2 2 
= “ + + Z + + + + 
A(3) = 290+326P+11P KK. 2K KZ KK 4 2K 2 K. 


2 2 
= + + + + + + + 
A(4) = 163+22P+3P 2K KZ KK 2 2K KZ KK, K, 
2 
=] +6P+P +2KK 42+K + 
A(5) 6 K, fs KK, 
A(6) = 1.0 
Taking the partial derivatives with respect to the system gain, 


AD(1) = z 


AD (2) =K2" +22 +2" 
t 
a AD(3) = KZ + 2K.Z + Kz 24 + 1.0 
AD(4) =K 2? + 2K Zt 2K EK. ET 6 
a a Le i 
AD(5) = 2K Z+K +K 
a a t 
AD (6) = 


AD(7) = 0.0 
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From the above definitions it is easy to imagine the complexity of the 
resultant equations when the bandwidth and sensitivity equations are 
expanded. These are practically impossible to solve without a computer. 
Using the computer program in Appendix D, a solution was obtained and 
is shown in the following pages. 

The solution of this system indicates that in order to meet the 
specifications the compensator zero of the double section of filter must 
be replaced on the right hand plane. By inspection of the compensated 
system roots the system is stable. This means that starting from a 
completely unstable system the result after compensation is a stable 
system satisfying all five specifications. A much more reasonable answer 
would have been if the compensator zero can also be moved to the left 
hand plane. This, however, is very difficult to meet with the given 


combination of compensator devices and required specifications. 
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INITIAL VALLES FCR UNKNCWN 


Ge SEQCCE C2 G.10G960 
Ge4G0Erb Co De 1UGlto 


AY 
go 


PARAMETERS ARE 


SPECIE ICATICAN FCR ROQT SENSITIVITY IS 


SPECIFICATICK FCR BANDWIDTH IS 


SPECIFICATICAN FCR ERRGR COEFFICIENT IS 


SPECIFICATICN FER CAMPING RATIC IS 


SPECIFICATICA FCR AATURAL FREQUENCY IS 


FCRWARD GATA IS 


TACHOMETER GAIN IS 


ACCELERATICA GAIN IS 


COMPENSATCR PCLE IS 


CCMPENSATCR Z7ERC IS 


NUMBER OF ITERATICAS IS 
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92375260 22 


92156840 21 


&.89858N-91 


9410010 09 


“92161190 1 


21 


~).2300UD 
ve 3000CD 


O11 


cyt? 
Tk 


Ue 200090 


a 
hae 


VL 


3.500000 34 


U.50NGID 2 


« 


TI.CALCULATEO SPECIFICATIONS 


ABS.VALUE CF CLTF AT RIS 9e79711D 06 


SENSITIVITY IS 92209000 0d 


ERRCR CCEFFICIENT IS 0e20000D O1 


ILTI/STABILITY INFCRMATICN 


COEFFICIENTS CF THE CeEe. IN ASCENDING ORDER 


3-146260 03 0.345910 63 9325740 03 


0.118780 03 046223D 52 0e719210 O1 


0.100000 01 ; 


THE CCMPENSATEC SYSTEM ROOTS ARE 


REAL PART IMAGINARY 


-0.71558D oc 0.421340 00 


-C.71558C0 GC 0.42134) GO 


-C.380450 ce “0.288780 01 


~C.38045C CC 0.28878D 01 


-0.250000 ¢1 “0.433010 01 


-0.2500CC O01 0.433010 01 
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Example 5.3 


Redesign the system of example 5.2 with a 3 db bandwidth at the 
frequency B = 0.34 and with a root sensitivity constant of 5 = 0.16. 
Inspection of the results shown in the following pages indicates that 
the given specifications are exactly met and the system is stable with 
all its poles and zeros on the left hand plane. Also it is worthwhile 
to observe the increase in system gain as the sensitivity constant 


was reduced. 
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INITIAL VALUES FOR UNKNOWN PARAMETERS ARE 


C.Ssagocn C3 0.20000D 06 
c.e70000D CC 0.50000D-01 


SPECIFICATICN FOR ROOT SENSITIVITY IS 


SPECIFICATION FOR BANDWIDTH IS 


SPECIFICATION FOR ERROR COEFFICIENT IS 


SPECIFICATION FOR DAMPING RATIO IS 


SPECIFICATION FOR NATURAL FREQUENCY IS 


TeUNKNCWN PARAMETERS 


FORWARD GAIN IS 
TACHOMETER GAIN IS 
ACCELERATION GAIN IS 
COMPENSATOR PCLE I§ 
COMPENSATOR ZERO IS 


NUMBER OF ITERATIONS IS 


74 


0.942910 03 


0-25084D 00 


0.47869D-01 


0-84389D 00 


0-66186D 00 


379 


9.50000D 06 
0-40000D 0C 


0.160000 00 


0.34000D 00 


0.20G00D 91 


0.50C00D 00 


0.500000 01 


ITeCALCULATED SPECIFICATIONS 


ARBS.VALUE OF CLTF AT B IS Ue7O07110 00 


SENSITIVITY IS 


0.16900D 90 


ERROR COEFFICIENT IS 92.2000CD O1 


TIT.STABILITY INFORMATION 


COEFFICIENTS CF 


0.619570 C3 


O.17557D C4 


g.10C00D 61 


THE COMPENSATED 


REAL PART 


-C€.74002D CO 


-C€.74002D C0 


-0.10485D Ol 


-C.2500CD G1 


-C.25000D Cl 


~0.42255D C2 


THE CoE. IN ASCENDING ORDER 


9.23763D 04 0-32005D 04 


0.358170 C3 02498230 02 


SYSTEM ROOTS ARE 


IMAGINARY 


-0.10591D 00 


0.165910 00 


0.0 


-0.43301D 61 


0.433010 Ol 


0.0 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The purpose of this paper was to use a digital computer to design 
a linear control system, which must satisfy up to five specifications, 
using algebraic methods. The feasibility of doing this was established 
and several examples showed that within certain limitations the method 
is applicable. 

Two computer programs were.written by the author using in both 
algebraic equations, developed by other authors in earlier work, to 
express specifications such as root locations, bandwidth, root sensitivity 
and error coefficient. The first program is accurate, fast and capable of 
handling root and error coefficient specifications to design a linear 
control system whose differential equations are available and can be 
expressed as a ratio of two polynomials in the complex:.frequency domain 
(transfer function). A fixed configuration of compensator devices con- 
sisting of a cascaded single filter and velocity (tachometer) feedback 
was used and the resultant nonlinear algebraic equations were solved 
by a very straightforward method of linearization. The program was 
written to handle uncompensated transfer functions up to nine order 
denominator and seven order numerator, however, implementation for 
higher order cases is apparent. On the other hand the program was 
written to give two solutions for a particular system if all the parameters 
are physically realizable. In addition, it was provided the possibility 
of getting only one solution or simply no solution when the parameters 


are not physically realizable. 


76 


As a consequence of increasing the number of specifications a 


different approach was needed for writing the second program. That is, 
the presence of specifications as root sensitivity and bandwidth made 

it impossible to use a linearization method to solve the resultant non- 
linear equations. Therefore, an interpolative procedure (Appendix C) 
was required and several subroutines and external functions sub-programs 
were written in conjunction with a rather general program based on the 
characteristic equation of the particular system to be designed. The 
examples showed the accuracy of the design and the dependence of the 
results on the set of initial values for the unknown parameters guessed 
by the designer. This set of initial values constitute the main dis- 
advantage of the technique (the method may diverge for poor initial 
guess for the unknown parameters). One way of implementing the initial 
guess could be by using some kind of search procedure as it is used in 
some professional subroutines (this is not shown in this paper). Another 
way used by the author, can be summarized as follows: 

a. Use program ROSDEP (roots and error specification design 
program) and find the unknown parameters required to meet root and 
error specifications using single section of filter in cascade and 
tachometer feedback as compensator devices. 

b. From the characteristic equation of the compensated system 
(actual problem) find the values of root sensitivity and frequency at 3 db 
bandwidth. This can be done using the same equations required as input 
for the actual program of design, FISDEP, by a very simple auxiliary program 


whose flowchart is indicated in Figure 6-1. 
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START 


Assume a 
value for 
‘ Ka 


ead values 
‘from ROSDEF 


Compute coef. of 
C.E. and partials 


ompute | 
Sensitivity | 


| 


Write Bandwidth 
Equations (BC) 


equal to 
0.707 + 


yes 


Write Bandwidth 
& Sensitivity 


Figure 6-1 


Flowchart showing sensitivity 
and Bandwitch calculations 
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om Compare the values obtained for sensitivity and bandwidth 


with the corresponding specifications and make reasonable adjustments 
according to the configuration of compensator devices required to use 
and some rules as: 

Ls The sign of the sensitivity constant is only a measure 
of direction. 

2. There is a tendency (verified for some systems) to 
increase the gain as the sensitivity decreases. 

The reason for the author to suggest this procedure is based on the 
fact that during the research he found that for several transfer functions 
studied the roots and error coefficient specifications were satisfied 
(even though the initial guesses were taken at random), but the sensitivity 
and bandwidth specifications were much more difficult to satisfy. There- 
fore, a Knowledge of what these quantities will be leads to a better 
guess and helps to insure the convergence of the procedure. On the other 
hand the total computational time using the suggested approach is 
increased only by about 30 seconds compared with the computational time 
required for a complete run taking directly the guessed values. 

In conclusion it was shown in this paper that a digital computer 
can be used to design linear control systems and that the method insures 
neither dominance of the given complex poles nor physical realizability 
of the system parameters, but it does insure that the specifications will 
be met and stability information can be obtained from the compensated 


system roots. Furthermore, the auxiliary program suggested to help in 
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finding the initial guess values to be used in program FISDEP (five 
specifications design program) provided a trivial method of computing 
root sensitivity and 3 db bandwidth for a root and error specifications 


design problem whose solution was obtained using program ROSDEP. 
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T8 


ba 


0.0 

0 .05C 
0.100 
0-150 
0-269 
9.25 
0.300 
0.350 
0-400 
0.450 
0.500 
9.550 
0-600 
92650 
0.700 
C.75f 
.80C 
9.850 
c.SCC 
9-950 
1.cCf 


U_LF) 
-1.609 
-1.0€2 
-1.0¢C0 
—-1.009 
-1.009 
1.000 
-1.000 
-1.009 
~1.0C0 
1.009 
-1.009 
-1.009 
~1.000 
-1.000 
-1.C¢9 
-1.006 
-1.009 
~1.009 
-1.069 
-1.009 
-1.609 


Ue) 
9.0 
0.0 
9.9 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
9.0 
9.0 
0.0 
0.6 
9.0 
9.0 
9.9 
0.0 
0.0 
0.9 


U,(F) 

-1.0990 
-1.9949 
—1.909 
-1.00f 
~1.995 
-1.090 
-1.9990 
~1.006 
-1.9096 
-1.090 
-1.999 
-1.000 
-1.099 
-1.999 
-1.900 
-1.0945 
-1.992 
~1.900 
-1.009 
-1.990 


“1.020 


APPENDIX A 


U,(F) 
-C.0 
9.2100 
0.200 
f.300 
.400 
0.5090 
0.600 
®.700 
C.806 
0.999 
1.009 
1.100 
1.299 
1.30¢ 
1.409 
1.560 
1.600 
1.799 
1.899 
1.994 
2-609 


TABLE OF CHEBYCHEV's FUNCTIONS (U, ( 


U,(F) 
1.000 
0.990 
32967 
9-910 
Je 849 
0.750 
9.2649 
9.519 
32360 
9.199 
9.900 
-9.219 
—7.440 
-3 699 
9.960 
-1.250 
-1. 569) 
-1.899 
2240 
-22610 
-3.900 


-1.88 
“1.071 
-1.900 
-0. 869 
22672 
0 6403 
-0.056 
0-375 
0.896 
1.513 
20232 
32059 
4.009 


~C 2546 
-$.313 
-9-059 
G.230 
0.2510 
0-774 
1.9900 
1.166 
1.246 
12214 
1.038 
0-688 
Q-126 
02682 
-1.778 
~3-262 
-52009 
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¥ 


0.0 

0.05¢ 
0-100 
0.159 
0.200 
9-250 
0.300 
0.2350 
0.460 
0-450 
9-500 
0.550 
0.600 
0.650 
0.76C 
0.750 
0.800 
0.85¢ 
0.900 
0.956 
1.000 


ULF) 
22000 
2.009 
2-000 
2.000 
2.000 
2.000 
2.000 
2-000 
2.000 
2.0C0 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2.000 
2-000 


CHEBYCHEV's FUNCTIONS (U; (¥)) 


ULF) 
0.0 

9.400 
0.800 
1.200 
1.600 
2-000 
22400 
2-800 
3-200 
3-600 
4.000 
42409 
4. 800 
5-200 
5-600 
6-000 
66400 
6-800 
7.200 
7.650 
8.000 


APPENDIX B 


FIRST DERIVATIVE OF SOME 


ULF) 
-4.900 
~3.940 
-3.760 
3.460 
~-32040 
2.509 
—1.840 
1.060 
-9.160 

0.86N 

2.000 


3.260 


4.2640 
62140 
7. 760 
92500 
11. 369 
13.340 
15-440 
17.660 
20.000 


UL(F) 
-0.0 
-1.192 
-~2.336 
-3. 384 
-4.288 
-5.900 
50472 
-5.656 
5.504 
4.968 
-4.000 
-2.552 
-0.576 
1.976 
5.152 
9.000 
13.568 
18.904 
252056 
320072 
40.000 


ULF) 
6.900 
5-761 
52956 
3.921 
22416 
92625 
-12344 
-3.359 
—-5 2264 
-6.879 
-8.000 
-8.399 
-7.2824 
~5.999 
22624 
20625 
10.096 
20.161 
33.216 
49.681 
79.009 


UF) 
9.0 
2.360 
4.484 
60149 
72163 
7.375 
6.693 
5.097 
20652 

0.2474 
4.090 
~7.514 
-10 460 
~12.125 
-11.621 
—7.875 
0.389 
14.663 
36.668 
68.371 
111.999 


UA(F) 
-8.000 
-T2406 
-5.695 
-3.076 
0.2121 
32469 
6e47TT 
8.641 
92494 
8.674 
6-900 
1.556 
—42212 
-10.391 
~152483 
-17-281 
12.735 
22290 
32.2715 
850315 
167.3998 


APPENDIX C 


The algorithm summarized in the flowchart of Figure C-l1 was used 
to write the subroutines to solve the nonlinear equations resultant from 
a five (or less) specifications design problem. "It is a variation of 
Newton's method, in which the partial derivatives are replaced by the 
multidimensional analogues of difference quotients. It produces super- 
linear convergence (of order 1/2 (1+ ket or approximately 1.62) to a 
simple zero, provided the given functions are twice continuously 
differentiable in a neighborhood of the zero, and the initial approximations 


are close enough." [6] 


K = Ktl 


Return 


+ 
Compute pk ; 


Figure C-l 


Flowchart showing the procedure 
used to solve a set of nonlinear equations 
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The capital letters indicated in the flowchart of Figure 3-1 refer to ‘4 


matrices (or vectors). To compute them the following definitions are 


required: 
ere thee x xi ak 
where 
lIx* ~ xk] || is the Euclidean norm of the vector X 
kik k_k kek k 
PP. PB grcrrecesecePyP Pr 
k, 2 kak kek k 
-(P,) PoP a: eee vee eePSP -P, 
k k,2 k,2 
H” = 0 ~(P)) + (P,) lee (e 
n-1 
0 | ae Sey" cer 
j n 


P =X -~X 
BX = £XK*) 
A=10"+p*) - BX 
eX =a7'p* 
pk - p*rk . 
Using these definitions andX =[x yzpr Mg and P = [dx dy dz dp ave A 


for a system of five equations and five unknown the following matrix 


can be written: 
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sdy/u sdxdz/vu sdxdp/wv dxdr/w -dx 


-sdx/u sdydz/vu sdydp/wv dydr/w -dy 
DS 0 -su/v sdzdp/wv dzdr/w -dz 
0 0 -vs/w dpdr/w -dp 
0 0 0 -w ~dr 
where 
2 2.1 
u = (dx + dy) ‘< 
2 2 
v= (ee + dy + dz V2 
WwW = Co + ay" + aoe + ap’) 1/2 
Z 
s = Gc + dy + ae" + an 5 ar?) 1/2 


The vector B and matrix A were computed using an external function (f) 
in which the nonlinear equations were set in terms of the unknown 
parameters. Likewise the vector E was calculated using an external 
subroutine for solving linear algebraic equations. 

Different subroutines were written for the cases of three, four and 
five simultaneous equations. In all of them a maximum of one thousand 
iterations were considered and indication of convergence was provided 
by using a counter which finally will be output of the main program. If 
the number of iterations in the output is indicated by a zero the 
procedure didn't converge and another set of initial values should be 


tried. Other parameters in the subroutines are self explanatory. 
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COEFFICIENTS OF THE COMPENSATED CHARACTERISTIC EQUATION FOR SOLe#2 
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